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1. Introduction 

Let (Af",g) be a smooth, compact n-dimensional Riemannian manifold with 
boundary and let A be the corresponding Laplace-Beltrami operator acting on 
functions. If the boundary is non-empty, we assume that either Dirichlet or Neu- 
mann conditions are imposed along 9M". 

Consider the operators x\ defined as projection onto the subspace spanned by 
the Dirichlet (or Neumann) eigenfunctions whose corresponding eigenvalues — 
satisfy Xj G [A — 1,A]. In the case that 9M" is empty, it was established in [10] 
that the following, best possible LP' — > estimates hold for x\'- 



(1-1) IIxaIIl^^l, 



< 



Recently, in [T] and [2], Burq, Gerard, and Tzvetkov established multilinear 
versions of these estimates, also under the assumption that the boundary of M is 
empty. To state these, suppose that A > /i > > 1. Then 

(1-2) \\x\f xMl^m) < A(^)||/||l2(m)||5IU2(a/) 

2ti — 3 

(1-3) \\X\f XtigXfhWL^iM) < (M'^)^~ll/llL2(Af)||g||L2(M)||^||L2(M) 

where in the first estimate 



A(m) 



With the exception of the logarithmic loss for n = 3, the linear estimate (|l.ip with 
g = 4 follows as a corollary of the bilinear estimate (|1.2p . by taking X — fi and 
/ = g. Similarly, the trilinear estimate p.3p implies (jl.ip with q — 6. Moreover, by 
taking h constant and = 1, (|1.3[) implies (|1.2p in case n = 2. For n > 4, however, 
the trilinear estimate can be improved by using p.2p together with the L°° bounds 
PTT|) on h. 

In the case where 9M" is nonempty, the issue of spectral cluster estimates is 
considerably more intricate. Here the Rayleigh whispering gallery modes provide 
examples of spectral clusters which concentrate in a A^^ x A ~ neighborhood 
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of a boundary geodesic (see Grieser [4]). These examples show that one cannot 
achieve linear spectral cluster estimates better than 



L2 




— ^ — 3n— 4 

fi^^ <q<oo 



(1-4) IIxaI 

3n-4 

The estimates ()1.4p were recently proven for dimension n = 2 in [9], along with 
partial results in higher dimensions. The question of whether or not they hold in 
general in higher dimensions remains open. 

In this work, we establish the following multilinear spectral cluster estimates on 
a general Riemannian manifold with boundary. We restrict attention to dimensions 
n — 2,3, since that is where our results are in some context sharp. 

Theorem 1.1. Let {M",g) and x\ be as above, with either Dirichlet or Neumann 
eigenfunctions, and let X > fi > i^. Then the following bilinear estimate holds 

(1-5) IIxa/XmSIU^CM) < A(/i) ||/|U2(M)||.g||L2(M) 

with A defined as 

/i3 n = 2 

IJ,i (log/i)5 n = 3 
In addition, the following trilinear estimate holds for n — 2,5 

(1-6) \\X\f Xf,9Xuh\\L2{M) < {l^l')~^\\f\\L^M)\\g\\L^M)\\h\\L^M) 

For n — 2, the bilinear and trilinear estimates imply the estimate (II. 4p . respec- 
tively for g = 4 and q = 6. Moreover, p.6p implies (|1.5p for Neumann conditions 
by taking h constant. For n = 3 this is no longer the case. However, the bilinear 
estimate for n = 3 implies (up to a logarithmic loss) the best possible L** linear 
estimate for manifolds with Lipschitz metric, by the examples of JS], and our proof 
in fact establishes Theorem 1 1.1 1 in this context. 



A(m) 



As in [^, a key step is to work on the double M of the manifold M, obtained 
by attaching two copies of M along the boundary, and taking coordinate patches 
along dM C M which agree with geodesic normal coordinates (x', a;„) on each copy 
of M. In these coordinates, the lift g of g to M is given by g^-'{x' , |a;„|), and hence g 
extends to M with a Lipschitz type singularity along dM. Dirichlet and Neumann 
eigenfunctions on M correspond to eigenfunctions on M which are, respectively, 
odd or even under a;„ —f —Xn- Hence, bounds on spectral clusters, and the 
multilinear analogs we consider, can be obtained by proving the same bounds on 
(M,g). 

The linear estimates of [9] were obtained by establishing mixed-norm L%Ll esti- 
mates on M for the evolution of a spectral cluster under the wave equation, and we 
follow a similar approach here. In that paper, the precise nature of the singularity 
of g along dM was used, and a microlocal decomposition of the cluster was made 
in terms of angle from tangent to dM . Estimates were obtained over small slabs, 
with size depending on the frequency and angle, and summing over slabs led to a 
frequency dependent loss for the estimates on M . 

In contrast, the results of this paper go through generally for the case of a 
boundary-free Riemmanian manifold with metric of Lipschitz regularity, as with 
the linear spectral cluster estimates of [7], or the Strichartz estimates of Tataru 
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[llj . As in those papers, we obtain estimates over small slabs with size depending 
on the frequency, and use a rescaling argument to reduce matters to obtaining 
estimates for metrics. We then use wave packet methods to obtain dispersive 
estimates, as in [5], [6] and [11], Summing over slabs then leads to a frequency 
dependent loss. 

2. MiCROLOCAL Reductions 

For the remainder of this paper, we assume M is a compact manifold without 
boundary, and g is a metric of Lipschitz regularity. The condition that /, h be 
spectrally localized can be relaxed, and we work instead with a quasimode condition. 
We state the condition for / here, the condition for g and h being analogous. For 
each local coordinate chart we write 

(2.1) gd^f + X^J^w, gd^^J2s'Hx)djdk. 

Given such an equation, we set 

lll/IIU - II/IIl^ + x-'\\dfh. + x-^d^fh2 + x-^wh. . 

If / = xa/, then for ip a smooth cutoff to a local coordinate system, the function 
ipf satisfies the equation ()2.1|) on M", and 

III0/IIU< ||/||l2(m). 

For the bilinear estimates it thus suffices to prove, for each coordinate chart, that 

II0/05||l^<A(a*)|||0/|||a|||</'5||U, 

and analogously for the trilinear version. 

By choosing appropriate coordinates, we may assume / satisfies (12. ip on M", 
with / supported in the unit ball, and 

||g^^ -<5*^ ||Lip(E") <Co, 

with Co a constant to be chosen suitably small. 

Let Sr = Sr{D) denote a smooth cutoff on the Fourier transform side to fre- 
quencies of size 1^1 < r. Let g^ — Sc^xg, for c to be chosen suitably small. Then 

\\{E-g>.)d'f\\L^<c-'X-^d'fU., 

and thus we may replace g by gA in (|2.ip at the expense of absorbing the above 
term into w, which does not change the size of |||/|||a- 

We next take a microfocal partition of unity, 1 = where each r(^) is a 

smooth symbol of order supported in a cone of small angle. By the Coifman-Meyer 
commutator theorem [3|, since g is Lipschitz 

[g,T{D)]d : L^{R") ^ L^{R"), 

hence T{D)f satisfies the equation with |||r(£))/|||A < |||/|||a • 

Since there are finitely many terms, we may replace / by T(D)f, which is no 
longer compactly supported, but is rapidly decreasing and smooth outside the unit 
ball. Without loss of generality we assume that r(^) is supported within a small 
angle of the axis. We similarly replace g and h by r'{D)g and r"{D)h, localized 
in frequency to small cones along general directions. 
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Letting x' — X2 in case of dimension n = 2, and x' = (x2, xa) in case of dimension 
71 = 3, we bound 

ll/5lU^ < ll/llLs=L^ hhiL-, 

XI a-' XI ^' 

Since r{D)f is rapidly decreasing outside the unit ball, it suffices to take the norms 
above over the ball of radius 2. Theorem If .H is then a result of the following. 

Theorem 2.1. Suppose that f satisfies the equation 

(2.2) g^dV + X^f = w 

Then the following hold, where the norms on the left side are over a bounded set 

(2.3) ll/IU^^L^°^ <Ai(logA)^|||/||U, n = 3 



(2.4) ll/IUjL^ 



< 



A^III/IIU, n = 2 
A^III/IIIa, n = 3 



Futhermore, if f{^) is supported in a small cone about the ^1 axis, then 
(2.5) ll/IL^^L^ <III/IIU 

■"^l x' 

Proof. We start by localizing / dyadically in frequency. Let 

/ = S,xf + (5,-u - Scx)f + (f - = f<x + fx + />A . 

Since [Scx,g\\d : — » L^, fx satisfies (12.21) . with |||i«a||U ^ lll/IIU, and similarly 
for />A and /<a. Furthermore, by the frequency localization of gA, each of w\, 
w<A, and w>A is also localized to the appropriate range of frequencies. 

A simple integration by parts argument (see the proof of Corollary 5 of 0) yields 
that, for c sufficiently small, 

A||/<aIU=(R") + IM/>aIU^(E")<|||/|||a. 

Elliptic regularity additionally gives the bound 

IMV>a||l^(M") <A|||/|||a. 

Sobolev embedding then yields each of the estimates (|2.3p - (|2.5p for /<a and 
Indeed, there is a gain of A 3 (log A) 2 in the estimate (j2.3p . and a gain of A 12 in the 
estimate (12. 4L for these terms. 



Consequently, we are reduced to establishing (|2.3p - (l2.5p for the term f\. We 
start with (12.51). Let V denote the vector field 



V = 2(ai/A) gA dfx + (AVI - (gA dfx,dfx))ri . 

Then 

div V = 2{difx) (divgA) • dfx + 2{d^x) wx - {{digx)dfx, dfx) 
Applying the divergence theorem on the set xi < r yields 

/ ^1 dx' < A2||/A||i.(R„) + |M/A||i.(M.) + ||ii'A||i2(M„) . 
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Since gA is pointwise close to the flat metric, we have pointwise that 

V^>^,\^^f,\^ + |X^\h\^-\^,,h\\ 
The frequency localization of f\ to |^'| < cA yields 

/ Vidx'>^f M/Ap + A^l/Apdx'. 

J xi—r ^ J xi—7' 

Consequently, 

^"^IM/a II L~L2, + II/a||l-l2^ < ||/a|U2(r„) + A"^||c?/a|U2(r„) + A"^||wa||l2(r„) 

< III/aIIIa, 

yielding estimate p.Sp . 

For the remainder, we assume that /(^) is localized to a small cone along the 
direction uj. In this case, the above argument yields uniform bounds over hy- 
perplanes of the form oj - x — r. In the proof of (|2.3p - ()2.4p we will use the following 
consequence. Suppose that Sn is a slab of the form lu ■ x £ I, where / is an interval 
of length |/| = R. Then 

(2.6) X-'\\dfx\\mSn} + II/a||l^(S„) < i?^lll/A|||A ■ 

We cover the bounded set on which the norms in (|2.3p and (|2.4p are taken by 
w A 3 slabs of the form uj ■ x E I, where |/| = i? = A" 3 . Then 

(2-7) II/a||lS,l- ^ sup ||/a||l5 r~(Sfl) • 

J- x' S Ji ^ 

We will establish the following result. Suppose that Qr is a cube of sidelength 
R = A^^, and Q}^ its double. Then 

(2-8) \\fx\\L'^^L^,{Qi,) 

<Cp(A)i?-^(||/A|U.(Q^)+A-l||#A||L^(QJ^)+i?A-imlU2(Qj^)) 

where 

C2(A) = A5(logA)5 , n = 3 c4X) ^ 

If we cover the slab Sr by disjoint cubes Qr, then by (12. 6|) we obtain 

II/aIIlSjL^csh) ^ Cp(A) |||/a|||a , 

and (1131) yields 

The estimate (12.81) is scale-invariant. Precisely, if we change x — > Rx, so that Q 
becomes a cube of size 1, and f\{R-) is frequency localized at scale RX = Aa , then, 
with /i = As, estimate (|2.8p is equivalent to the following 

(2-9) ||/mIIls^l-(q) < Cp(Ai)( ||/mI|l2(q*) +^"i||d^||i2(Q*) +/i"lu'Mlli^(Q*)) ■ 
Here, /^(a;) = f\{X^^x), which satisfies the equation 

with gp(a;) = g\{X^^x). Observe that 

(2.10) ||rfgMl|z,~ < coA"^ = cofi^i , 




6 



MATTHEW D. BLAIR, HART F. SMITH, AND CHRISTOPHER D. SOGGE 



hence 

with g^i/2 = S'f,2^i/2g^ . Thus satisfies the equation 

(2.11) g^^/2d% + ^l^f,,^w^, 

with the right-hand side of (|2.9p of comparable size. 

It follows from ()2.10p that the metric g^i/2 is of regularity C^. Indeed, 

\\g';.,,_-p'\\c^<co. 

To establish (|2.9p . we may thus use techniques developed to establish dispersive 
estimates for operators of principal type with coefficients. We follow below 
the path through squarefunction estimates for solutions to a first order hyperbolic 
equation, by introducing a time variable, as in [6]. It should also be possible to 
establish the dispersive estimates directly for (|2.1ip using methods of [5]. 

Let 
Then 

\\p{x, Dff^ + g^i/2dVMllL2(R") < M II/mIU^ ■ 

Thus, 

(m + P{x, D)) - p{x, D)) ff^^Wf,, 
with the error harmlessly absorbed into w^. The operator /i + p(x, D) is elliptic on 
the frequency support of (/i — p{x^ D)^ f^^ hence wc may write 

{p-p(x,D))ff, = p^'^Wf, , 

with 1 11/11 |a still of comparable size. Finally, let 

u{t, x) = e-^*^^(x) , F = fi-^e-^'^w^ . 

Then 

{dt+ip{x,D))u^ F , 

and it suffices to show that 

(2.12) ||u||l£jL~L2(qx[04]) < Cp(/i)( ||u||l~L2([o,i]xR") + ||^'||l1L2([04]xR")) . 

Our proof of (|2.12p follows very closely the proof of the linear spectral cluster 
estimates in [B]; we outline just the main steps here. Following [51 §3], consider the 
wave packet transform of u^, 

(T^u) (t, X, e) = J e-*<«'"-"> (z - x)) u{t, z) dz , 

where is a real, even Schwartz function, with ||0||l2 = {2tt)^^ ^ and with Fourier 
transform supported in the unit ball {|^| < 1} . Then 

dtT^uit, = {d^p{x, 0-dx- dxp{x, ^ ■ d^ Tf,u{t, x, C) -f G{t, x, , 

where G{t,x,S,) = unless ^p < |^| < 2p, and 
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Let xt denote the canonical transform on M^"^ — r*(M") generated by the Hamil- 
tonian flow of p. Thus, Xtix,^,) — "f{t), where 7 is the integral curve with 7(0) = 
(x,^). Then we have 



Thus, T^u{t,x,S^) is an integrable superposition over r of \t>r multiplied by a 
function invariant under the Hamiltonian flow of p. 

Since u{t,x) — T*(T^uj{t,x,£,), it sufRces to show 

(2.13) II WIIl£,l~l?(r..x[o,i]) < Cp(a*)||/||l^,^ , 
where 

{Wf){t,x)=T;{fox-t){x). 
This is in turn equivalent to the following bounds 

(2.14) II W*F||^.^,.^.(Q,[o,l]) < ll^^ll<...,^.(gx[0,l]) • 

The operator WW* has an integral kernel K which is highly localized to a 
neighborhood of the light cone, with the dispersive rate of decay away from the 
origin, see O (3.11)], 

(2.15) \K{s, y;t,z)\<^l''{l + ^l\y,-z,\)''^ + 0) - - i| | ) , 

with d{y, z) the distance of y to z determined by p. 

We remark that in [6j this estimate was established assuming the kernel was 
microlocalized near the ^1 axis. That assumption, however, was necessary for L^- 
energy estimates, not the above dispersive estimates. Indeed, the proof of [H (3.11)] 
establishes (|2.15p with \yi — zi\ replaced by \y — z\, since |i — s| ~ |y — ^1 on the 
light cone, and hence holds without any assumption of conic micro-localization. 

Estimate (|2.15p implies that, for each {yi,zi), 



K{s,y;t,z)v{t,z')dtdz' < /i""^(l + fi\yi~ zi\ 

y ^ 

For n = 3 and p = 2, estimate (|2.14l) follows from 



WW Tl T.2 . 



L\Lf 



/ Ai^(l + M|yi " ^i| ) ^<AilogM. 

J\z-, l<2 



l\zi\<2 

For p = 4, (|2.14p follows from the Hardy-Litlewood-Sobolev inequality, together 
with the bound 

M"~i(l + /i|yi-zi|)-'^ <M"-^|t/i-zir^. □ 
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